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Abstract 

Free energy as a function of temperature and the bubble radius is determined 
for spherical bubbles created in cosmological first order phase transitions. The 
phase transition is assumed to be driven by an order parameter (e.g. a Higgs 
field) with quartic potential. The definition of the bubble radius and the cor- 
responding generalized, curvature-dependent surface tensions are discussed. In 
the free energy expansion in powers of the inverse radius, the coefficients of the 
curvature term and the constant term are also calculated. 
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Recently the detailed mechanism of the electroweak phase transition, which 
depends on the free energy of the nucleating bubbles, has gained much interest [|], |], [3j , 
mainly because it is believed that this transition has had a significant effect on the 
baryon asymmetry of the Universe. The purpose of this Letter is to study the free 
energy of bubbles which are created in cosmological first order phase transitions. In 
particular, we study its expansion in 1/-R, the definition of R (the bubble radius) and an 
independent surface tension. The results can also be used to discuss the observation |4j] 
that in the MIT bag model small hadron bubbles exist already above the quark-hadron 
phase transition temperature. 

The phase transition of the electroweak theory, and of the possible grand unified 
theory as well, is generally thought to be driven by a Higgs field. The theory could 
contain several Higgs fields, but one combination can be used as an order parameter. 
For QCD, the choice of the order parameter is not so clear, since the theory has no 
classical potential. For example, the energy density could act as an effective order 
parameter field ||. It should be noted here, that even though we consider only first 
order phase transitions in this Letter, in reality the order of the electroweak or the 
quark-hadron phase transition is not known. 

We begin by writing the action for the bosonic order parameter field in the high 
temperature approximation: 



S = ^ MV 



i(V0) 2 + V(0,T) 



The effective potential, which is commonly used to describe a first order electroweak 
phase transition is || [l], ||, |[ 

T) = i 7 (T 2 - T 2 )0 2 - iaT0 3 + . (2) 

The potential for the energy density, which was used in || for QCD, is also quartic. 
This potential V(<f), T) in eq. (0) should be regarded as a phenomenological expansion 
for the effective order parameter, valid in the vicinity of T c . 

First, we will describe the relevant thermodynamical properties of the potential 
in eq. (0), based on 0. In the thermodynamical limit, the transition occurs when the 
heights of the two minima are equal. This happens at the temperature T c , 

T c = ; T ° s • (3) 

In cosmology, the transition takes place after some supercooling at a temperature lower 
than T c . It is convenient to express the temperature dependence of several quantities 
by using the function A(T), 
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which satisfies A (To) = 0, A(T C ) = 1. 

The pressure in the low temperature (broken symmetry) phase minus the pres- 
sure in the high temperature (symmetric) phase is 

ApCO = ^T 4 {^A(T) 2 - ~A(T) + 1 + [1 - 8 -X(T)f^ . (5) 

This pressure difference is positive for temperatures smaller than T c . At T c , the corre- 
lation lengths in both phases are equal to 

(6) 

The action in eq. ([I]) determines two additional thermodynamical quantities which are 
useful for us; namely the latent heat 



4cT7 



T % 2 (7) 
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and the thermodynamical or planar surface tension 

The latter quantity is denoted by cr p , since later on we will define a generalized surface 
tension a(R). 

The latent heat coming from the decrease of the effective relativistic degrees 
of freedom is included in eq. (|^). The constants T , 7, a and A are to be chosen so 
that the potential quantitatively correctly describes the phase transition. These four 
constants can be expressed in terms of the more physical quantities T c , l c , L and <r p 0: 

rp _ T c _ L+6a p /l c 



a 



P'C-'c 

(9) 



Values of Tq/T c which are smaller than a/5/3 ~ 0.75 are not very natural 0. 
However, when the potential is required to be valid only in the vicinity of T c , also 
smaller values are possible. Then, if one requires that the vacuum expectation value 
of the low temperature phase increases with decreasing temperature near T c , it follows 
that T /T c > l/y/3 w 0.58 . 

The decay rate of the metastable vacuum can be calculated using standard 
methods || |9|, |6|]. Probability of tunneling per unit time per unit volume in the high 
temperature approximation is 

P{T) = a (T) e~ s " {T) , (10) 
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where the pre-exponential factor a(T) is expected to be of the order of T c 4 . The critical 
action S cr (T) is the value of the action ([!]), with the appropriate boundary conditions, 
in the 0(3)-symmetric extremum. More physically speaking, it is the free energy 
F CT (T), which is needed to form a critical bubble, divided by the temperature. 

The critical action can be analytically approximated in the small and large 
relative supercooling regimes ||, and in the full range it has been calculated numeri- 
cally [fj, |3|. In terms of the function A(T) from eq. (f|), 

ci[ >~ T ~ 3 5 7r A 3 / 2 (1 - A) 2 ' 1 } 

A fit to the function /(A) with an —0.1 . . . +1.6% error was presented in |3[]. Here we 
give for /(A) a four-parameter fit with an accuracy of ±0.1% : 

/(A) = A 3/2 (oq + a x \ + a 2 A 2 + a 3 A 3 + a 4 A 4 ) ; (12) 

a = 15.63628, a x = -18.03398, a 2 = 2.39731, 
a 3 = -0.86504, a 4 = 1.86543. 

This fit is exact in the small relative supercooling (A = 1) limit. The function /(A) 
is shown in fig. [I], together with /(A)/A 3 / 2 . It is interesting to note that the latter 
one, which is relevant in the large relative supercooling (A = 0) limit 0, is an almost 
straight line over the full range. 

Now we will turn to the main issue of this Letter: study of the Helmholtz 
free energy of a spherical bubble of the low temperature phase nucleating in the high 
temperature phase, and a discussion of what effects differing definitions of the bubble 
radius have. In the theory of fluid surfaces, the dependence of the surface tension on the 
definition of the bubble radius is a subtle issue |TDfl . With the exact calculation of the 
extremal configuration from eq. (|l|) the effect of changing the definition can precisely be 
evaluated. One should bear in mind, however, that when discussing bubble formation 
only the free energy has real physical importance; surface tension as well as curvature 
coefficient are only auxiliary quantities. 

Expanded in powers of 1/R, i.e. in the large R limit, the free energy can be 
written as 

4:71 

F(R, T) = Ap(T)R 3 + 47ia f (T)R 2 + 8irj f (T)R + 16n6 f (T) + ... , (13) 

3 

where R is the bubble radius according to some specific definition. Note that this 
general bubble free energy is a function of both R and T, whereas the free energy of a 
critical bubble, F CT (T), depends on one variable only. They are related by 

F cr (T) = F(R cr (T),T), (14) 
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where R cr (T) is radius of the critical bubble. 

It is expected that near T c the relative change of the functions 07 (T), Jf(T), 
Sf(T) is slow, in contrast to the pressure difference Ap(T), since properties of the 
interface do not change dramatically with temperature. In cosmology, where we usually 
are discussing bubble formation only near T c , we can treat these functions as constants. 
This assumption enables us to solve the general free energy F(R,T) from the critical 
free energy F CT (T) and the pressure difference Ap(T). Now the free energy is for large 
R, in the vicinity of T c 



An 

F(R,T) = Ap{T)R 3 + 47ra p R 2 + S^R + 16n5 + . . . 
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Am 

= -—Ap{T)R 3 + Ana{R)R 2 , 



(15) 
(16) 



where the latter form defines the generalized surface tension (see also [fTI]) 



a(R) 



AttR 2 



Att 

F(R,T) + -jAp(T)R 3 



(17) 



The critical free energy F CI (T) determines a(R cr (T)), assuming that Ap(T) is known. 
If R cr (T) is known as well, we can solve for a(R). And from a(R), the coefficients in 
eq. (|15D can be solved: 



(Jr, 
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a{R)U = 
1 da(R) 



' 



2 d(J 
1 1 d 2 a{R) 



(18) 
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We employ two different definitions for the bubble radius. The first one is 
Laplace radius R^ defined by using Laplace's relation 

2a L {R Lcr {T)) 



RLcr(T) 



(19) 



Ap(T) 

The second one is the tension radius Rt, defined as the distance at which the gradi- 
ent density term of the critical action reaches its maximum. The definition of Rt is 
illustrated in fig. ||. 

From Laplace's relation, it follows that the radius of the critical bubble and the 
generalized surface tension in terms of the critical free energy are 



Rlct{T) 

OL(i?Lcr(T)) 



3 F(R hcI {T),T) 



2tt 
3 

167T 



Ap(T) 
F(R Lcr (T),T) Ap(T) 2 



(20) 
(21) 
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In the maximum tension definition, the radius of the critical bubble, Rt C t(T), is ob- 
tained from the numerical solution of the extremum differential equation for the action 
in eq. (|l|) for different values of A. The curvature-dependent surface tension cttC-Rtct^)) 
is then calculated directly from eq. (fPTD. 

The two definitions for the bubble radius are compared in fig. [3]. One notes 
that always i?Tcr > Rhcv In the small relative supercooling limit, the two definitions 
coincide. On the other hand, in the large relative supercooling limit they differ even 
qualitatively: i?Lcr vanishes, whereas i?Tcr goes to infinity. However, this limit is beyond 
the validity of the expression (15,16) for the free energy, in which all the temperature 



dependence comes from the pressure difference. 

Both of the generalized surface tensions, <7l(.R) and ctt(-R), are plotted versus 
the inverse radius l c /R in fig. [|. For large values of the bubble radius the two surface 
tensions coincide. In this limit, depending on the value of the constant T /T c , the 
surface tensions can either decrease or increase with R, implying that the curvature 
coefficient can be positive or negative, respectively (see fig. |5|). When the bubble radius 
gets smaller, the two surface tensions differ increasingly. The curves for <Jt(R) turn 
backwards, because with decreasing A, i?Tcr when expressed in physical units possesses 
a minimum, after which it begins to increase (see fig. ^). For small values of the radius, 
ctl(-R) has to decrease with R, in order to compensate the faster decrease of the term 
Ap(T)R 3 in the free energy expansion in eq. (0). And vice versa, ot(-R) must increase 
with decreasing A in the large relative supercooling regime. The two-valuedness of 
o"t(-R) does not imply that the distance of maximum tension would be a bad definition 
for the bubble radius. The radius R^ agrees well with the intuitive picture of the 
bubble, and furthermore, as will be discussed later, its minimum value can be used as 
an estimate for the break-down distance of the free energy expansion. 

The coefficients in the free energy expansion in eq. flT5] ) were determined from 
eq. (|TH), employing both definitions for the bubble radius. Within numerical accuracy, 
the results were equal. The reason for this is that both of the radii, and both of 
the generalized surface tensions as well, coincide in the limit of infinite bubble radius. 
Below, results for the coefficients are presented only for the case where the radius is 
defined by Laplace's relation. 

We write a^(R^ CT (T)) and IJ R^ cr (T) as Laurent expansions of A around A = l, 
which corresponds to infinite radius. We find that the coefficient of the curvature term 
is 
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where 



2 f(l) u 
9 9 2 



a p l c w (l .49 ± 0.01 - ^\ a p l c , (22) 



9A 7 -2a 2 2 I T 2 I v ' 
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In eq. (|22"D , f'(l) denotes d/(A)/dA at A = 1, and its numerical value is computed from 
eq. (|12"D . For the coefficient of the constant term in eq. (|i~5"D we find 



The coefficients 7 and 5 are shown in fig. [5[ If the temperature T is near T c , the 
cofficient 7 is in units of a p l c close to unity, and S almost vanishes. If T /T c is in the 
vicinity of 0.58, the zero-point of 7, the constant term dominates over the curvature 
term in the free energy expansion even for rather large bubbles. When Tq/T c decreases 
from ~ 0.4, 7 decreases and 5 increases rapidly. 

The magnitudes of the different terms in the truncated free energy expansion are 
compared in fig. || From the figure we can conclude that if T /T c > 0.5, the truncated 
expansion for the free energy of a spherical bubble in eq. flT5] ) breaks down when the 
bubble radius is of the order of 2/ c . If T /T c is smaller, the truncated free energy 
expansion is valid only for much larger bubbles. The distance where the truncated 
free energy expansion breaks down is for different values of Tq/T c roughly equal to the 
smallest tension radius of the critical bubbles. 

Finally, we will investigate if the low temperature phase bubbles can exist at 
temperatures higher than T c . Here the full expression for the free energy, eq. (|HJ), is 
used instead of the truncated one. Now the volume term in this equation is positive. 
Therefore the necessary condition for the existence of stable low temperature phase 
bubbles is that the surface term (i.e. the latter term) not be a monotonically increasing 
function of the radius R. The surface term is plotted in fig. [7] for two values of T /T c , 
and it behaves qualitatively in the same manner for all other values of Tq/T c . In 
spite of the two-valuedness of ctt(-R) it is clear that the behaviour of the surface term 
4:7t<j(R)R 2 shows no indication of the existence of low temperature bubbles above T c . 

- The expression for the free energy in eq. (|16D should not, however, be used for too 
small values of the bubble radius, because qualitatively it does not describe critical 
bubbles of large relative supercooling. The qualitative description breaks down when 
the bubble radius is again of the order of the smallest value of -Rtct- 

The result that the magnitude of the generalized surface term grows with the 
bubble size is very natural (even though u{R) is a decreasing function of the bubble 
radius for large bubbles if T /T c > 0.58, the combination <j(R)R 2 is an increasing one). 
However, in the MIT bag model calculation by Mardor and Svetitsky Q it was found 
that the free energy as a function of the bubble radius has a minimum for T > T c . 
Clearly the conclusion must be that either the spherical MIT bag model, or any order 
parameter model with a potential as the one in eq. (|2]), describes the quark-hadron 
phase transition qualitatively incorrectly. 
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I {/"(I) - 2/'(l) - 6 + 2u [5/(1) - 1] + 27u 2 } a v l\ 



(24) 




(25) 
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To summarize, we have discussed bubble free energy in cosmo logical first order 
phase transitions which can be described by the quartic potential. Numerical function 
for the free energy of critical bubbles with an ±0.1% accuracy has been presented. 
Two different definitions for the bubble radius have been studied: Laplace's relation 
and the distance of maximum tension. The corresponding curvature-dependent surface 
tensions <xl(-R) and ctt(R) have been calculated, and expanded to second order in 1/R 
around infinite bubble radius. It has been shown that for typical values of the constant 
T /T c a(R) is for large R a decreasing function (positive 7), but that in principle it 
can also be an increasing function (negative 7). No indications have been found of 
the existence of low temperature phase bubbles above T c . The assumption that all the 
temperature dependence in the free energy comes from the pressure difference between 
the two phases made it possible to solve the general free energy F(R, T) of spherical 
bubbles from knowledge of the critical bubbles only. 
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Figure 1: Non- divergent part of the critical action. Solid curve is /(A), as defined in 
eq. ([TT]) , and dashed curve is [/(A)/A 3 / 2 ]/6 . 

Figure 2: Definition of R^. Solid curve is the gradient part of the action density of a 
critical (A = 0.9)-bubble and dashed curve is the potential part, the first and the latter 
term in eq. ([!]), respectively. The dimensionless radius is R' = MR, where M stands 
for the bosonic mass in the potential, M 2 = 7(T 2 — TJ 2 ). The free energy density is 
shown in units of TM 3 a/\ 3 ^ 2 (note that since A is fixed, so is T/T c ). From the figure, 
we can read that for A = 0.9, i?^ cr ~ 7.3 . 

Figure 3: Radius of the critical bubble for Tq/T c = 0.58, according to the two defini- 
tions. Dashed curve is R^cr/h, and solid curve Rt ct /I c . 

Figure 4: The generalized surface tension from both definitions, o"l(-R) and ctt(-R), 
plotted versus the inverse radius l c /R. Dashed curves show ctl(-R)/ct p , and solid curves 
a^(R)/a v . Upper curves are for T /T c = 0.99, and lower curves for T /T c = 0.58 . 

Figure 5: The coefficients 7 and 5 in the free energy expansion in eq. (|T5|), for different 
values of T /T c . Dashed curve is 7/(<7 p Z c ), coefficient of the curvature term, and solid 
curve is 5/(cr p Z 2 ), coefficient of the constant term. The width of the curve for 5 shows 
the uncertainty in eq. (p5|), due to numerical determination of f"(l) (and to much less 
extent, of /' (1)). 

Figure 6: The distances at which magnitudes of the curvature term and the constant 
term equal the surface term in the truncated free energy expansion of eq. ( |i~5"D , plotted 
versus T /T c . Dashed curve shows the distance R 1 jl c at which 4na p R 2 = 8tt\j\R^, and 
solid curve the distance Rs/l c at which 47ia p Rj = 16ii\5\. 

Figure 7: Logarithm of the surface term in the free energy expansion in eq. (|T6D, 
log[47rcr(i?)-R 2 /(cr p / 2 )], as a function of the bubble radius R/l c - In the same manner as 
in fig. |], dashed curves represent 4it<tl(R)R 2 , solid curves A^a^{R)R 2 ] upper curves 
are for T /T c = 0.99, lower for T /T c = 0.58 . 
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